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Sigma meson in pole-dominated QCD sum rules 
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The properties of <r(600) meson are studied using the QCD sum rules (QSR) for the tetraquark 
operators. In the SU(3) chiral limit, we investigate separately SU(3) singlet and octet tetraquark 
states as constituents of the a meson, and discuss their roles for the classification of the light scalar 
nonets, <r, /o,ao 5 an d k, as candidates of tetraquark states. All our analyses are performed in the 
the suitable Borel window which is indispensable to avoid the pseudo peak artifacts outside of the 
Borel window. The acceptably wide Borel window originates after preparing the favorable set up of a 
linear combination of operators and the inclusion of the dimension 12 terms in the OPE. Taking into 
account for the possible large width, we estimate masses for singlet and octet states as 700 ~ 850 
MeV, 600 ~ 750 MeV, respectively, although octet states have smaller overlap with the pole than 
singlet state and may be strongly affected by low energy scattering states. This splitting of singlet 
and octet states emerges from the number of the qq annihilation diagrams, which include both color 
singlet annihilation processes, qqqq — > (qq)i and color octet annihilation processes, qqqq — > G(qq)s- 
The mass evaluation for the a meson gives the value around 600 ~ 800 MeV which is much smaller 
than the mass obtained by 2-quark correlators, 1.0 ~ 1.2 GeV. This indicates a state has the large 
overlap with the tetraquark states. 

PACS numbers: 12.39.Mk,11.55.Hx,11.30.Rd 



I. INTRODUCTION 



The structure of scalar mesons is a long-standing prob- 
lem in hadron spectroscopy In contrast to the other 
hadrons, two flavor nonets appear around 1 GeV in the 
scalar meson spectra. The systematic classification of 
scalar mesons sheds light on the nonperturbative aspects 
of QCD. Especially, the lighter scalar nonet, the isoscalar 
cr(600), / (980), the isodoublet k(800), and the isovec- 
tor ao(980), arc candidates for exotic hadrons (exotics) 
which have more complicated components than simple 
mcsonic qq or baryonic qqq structure. The naive qq 
assignment of the constituent quark picture for these 
mesons, a ~ {uu + dd), fo ~ ss, k ~ us, clq ~ (uu — dd), 
may lead the unrealistic mass ordering, m(cr) ~ m(ao) < 
m(n) < m(/o), and suggest heavier masses due to the P- 
wave angular excitation between two quarks, which typ- 
ically costs 500 MeV, for the J p = + quantum number 
of the scalar mesons. It is believed that these qq assign- 
ments are most likely realized in the nonet above 1 GeV 
rather than the nonet below 1 GeV. 

The masses of the light scalar mesons are explained 
by several pictures. One of them is the four quark pic- 
ture proposed by Jaffe Q. In this picture, assuming 
the quark contents with the ideal mixing of the flavor as 
a ~ (ud)(ud), fo ~ [(ds)(ds) + (su)(su)}, n ~ (ud)(su), 
ao ~ [(ds)(ds) — (su)(su)], one can naturally explain the 
observed mass ordering. For a meson, quenched lattice 
calculations also support this assignment @, ||. Quali- 
tatively, one of the key ideas for the considerable mass 
reduction of the light scalar mesons below 1 GeV is to 
employ the possible strong diquark correlation originated 
from the chromo- magnetic interaction 0] • This brings us 
to an interesting possibility that the scalar mesons are 



considered as the good place to investigate the strength 
of diquark correlation, which may provide a useful build- 
ing block to understand the hadron spectra [f| and for 
the further applications to the hot/dense quark matter 

Mixing of the two and four quark components is also 
employed to explain the scalar meson spectra. For ex- 
ample, the scalar nonet above 1 GeV have larger mass 
than one expected from the qq picture in the conven- 
tional quark models, as seen in the anomalous spin-orbit 
splitting, ai(1230) < a 2 (1320) < a (1450), in contrast to 
the charmonium mass splitting x c o(3414) < x c i(3510) < 
Xc2(3556). This mass ordering could be explained by a 
mixing scenario of the two and four quark components, 
in which, as a result of the level repulsion, the mass of 
the four-quark dominate state gets reduced, while the 
two-quark dominated state is pushed up 0] • 

For the isoscalar sector, a and fo, Narison discussed 
another possibility [|[ invoking the QCD sum rule (QSR) 
[9l.[Tcl| and some low energy theorems [ll|. In his scenario, 
the glucball degrees of freedom come into play and the 
strong qq — G 2 mixing leads the considerable mass re- 
duction below 1 GeV although unmixed glueball and qq 
are expected to be relatively heavy, ~ 1.6 GeV [l2j and 
~ 1.4 GeV [U, respectively, in the quenched lattice cal- 
culations. The strength of the qq — G 2 mixing has been 
studied in various approaches, but the strong mixing has 
not been confirmed yet [l3T | . 

The hadronic picture for the scalar meson were also 
widely discussed. The a, «(800), /o(980) and a (980) are 
dynamically generated as quasi-bound resonance states 
in scattering of the Nambu-Goldstone bosons based on 
chiral dynamics with an appropriate treatment for restor- 
ing unitarity in the scattering amplitudes [3] . 

The a meson itself is an attractive subject of contem- 
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porary nuclear physic. The a field is the origin of the 
attractive part of the nuclear force in the intermediate 
energy region flEl ] and can be a chiral partner of the 
pion being a possible soft mode in the chiral restora- 
tion Its existence of the a meson had been a long- 
standing problem, but recent experemental results such 
as D + — > 7r + 7r + 7r~ 17], J/ip — > ujtt + tt~ 03, and disper- 
sion analysis employing Roy equation [l9( support the 
existence of the a pole with mass 440 ~ 540 MeV and 
width 250 ~ 540 MeV. 

Yet the structure of the scalar mesons is not conclusive. 
We have several pictures of the constituents of the scalar 
meson with their admixtures. One of the purposes in this 
paper is to discuss the relevant constituents and the mix- 
ing scheme. They can be investigated through compari- 
son of the several correlators utilizing in the QSR 0, [l(| , 
which relates the nonperturbative aspects of QCD to the 
hadronic properties. In this work, our main consideration 
is on the four-quark (4q) picture for the light scalar nonet 
and we employ the tetraquark operators to obtain large 
overlap with the tetraquark states. We will show that the 
correlators of the tetraquark operators include not only 
four-quark connected diagrams but also disconnected di- 
agrams, such as the (qq)i and (qq)gG in the intermedi- 
ate states. These mixing effects can be studied from the 
difference between SU(3) flavor singlet and octet corre- 
lators. In fact, the number of the annihilation diagrams 
is four times larger in the singlet case than in the octet 
case. In our operator case, this leads larger low energy 
enhancement in the flavor singlet case. Throughout this 
paper, to avoid the complication from the current quark 
mass effects, we take the SU(3) chiral limit, in which the 
singlet and octet states completely decouple. 

For the tetraquark analyses with QSR of the a me- 
son, we have several technical issues to be solved. The 
most important point is that the resonance mass should 
be extracted within so-called Borel window in which the 
QSR works with better accuracy. The Borel mass M is 
introduced to the QSR by a derivative operation Lm to 
the dispersion relation of the correlation function II as 

L M n°P e {-Q 2 )= / ds e-^-MI^s), (1) 
Jo n 

where Lm = d n /(dQ 2 ) n taking n — ► oo and Q 2 — > oo 
being fixed by Q 2 /n = M 2 . The correlation function 
jjope - m ^ ne j e fl. nanc [ gi^g j s calculated by the opera- 
tor product expansion (OPE) and IP in the right hand 
side is expressed by hadronic contributions. The tower of 
successes in reproducing the meson and baryon spectra 
is attributed to appropriate and careful applications of 
QSR in suitable Borel mass region, i.e., Borel window. 
Only within the Borel window, one is allowed to extract 
the low energy properties of the spectral integral. 

The main technical problem for the multi-quark QSR 
is that the setting up the Borel window becomes quite dif- 
ficult. In contrast to the usual meson and baryon cases 
[20| , since correlation functions of multi-quark interpolat- 
ing fields show slow convergence of OPE and unwelcome 



high energy contamination dominates the spectral inte- 
grals. As we will emphasize in Sec lII Al and III BL this 
leads difficulties to extract low energy properties of the 
spectral function and we are often stuck with the pseudo 
peak artifacts outside of the Borel window. 

All these issues can be solved by use of suitable in- 
terpolating fields and inclusion of the OPE terms up to 
dimension 12 (dim. 12) for the tetraquark operators. Af- 
ter that, the Borel windows are found and we can in- 
vestigate the physical quantities within the Borel win- 
dows. To complete solving these issues, we would like to 
add one more analysis for the a meson in QSR with the 
tetraquark operator despite tetraquark operator analysis 
for the a meson, although there have been many works 
done in the past [2ll.[22l|. 

Another important aspect when we discuss the a me- 
son is its possible large width and the effects on the mass 
evaluation of a. In Sec lII C[ wc discuss the Borel trans- 
formed Brcit-Wigncr type spectral function and test how 
they behave in the effective mass plot. We find the sta- 
blity is moderate even up to ~ 400 MeV width. As 
related to this, we re-examine the usual criteria to fix 
Sth- They are reflected to the estimations for the physi- 
cal quantities of the a meson. 

The organization of this paper is as follows. In Sec. II, 
we explain the basic concepts of QSR and illustrate the 
typical pseudo-peak artifacts which we often encounter. 
The importance to set the Borel window is emphasized. 
The effects of the width on the effective mass plot are 
illustrated, and the criterion to fix Sth is re-examined. In 
Sec. Ill, we discuss the flavor singlet and octet states in 
the chiral limit. We also argue that both color singlet 
and octet annihilation diagrams are source to split the 
singlet and octet states. This splitting may be impor- 
tant to classify fo, K, oq even after taking realistic quark 
mass, which will be discussed in the subsequent paper. In 
Sec. IV, we show the results of the Borel analyses for the 
sigma meson together with the singlet and octet states. 
At first we introduce the tetraquark operator used in the 
analyses as a linear combination of two types of diquark- 
diquark local operators. We discuss the criteria to deter- 
mine the mixing angle of the operators. After fixing it, 
we perform the Borel analysis. For the singlet state we 
obtain the mass value 700 ~ 850 MeV, while, for octet 
state, we estimate the mass 600 ~ 750 MeV, although 
the residue of the octet state is much smaller than the 
singlet state and thus may be affected by low energy scat- 
tering background. Finally we examine the a meson as 
superposition of singlet and octet states and estimate its 
mass 600 ~ 800 MeV within error of width effects. The 
Sec.V and VI are devoted to discussion and summary, 
respectively All the details about the OPE terms and 
operator dependence are summarized in the Appendix. 
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II. THE BASIC CONCEPTS OF QSR AND 
POSSIBLE PSEUDO-PEAK ARTIFACTS 

In this section, we start with a brief review of the basic 
concepts of QSR, especially emphasizing the importance 
of the Borel window. Definitions of the terminologies and 
notations used in the later analyses are given in Scc lII Al 
In Sec lIIBl we illustrate a typical artifact in the sum 
rules, pseudo peak artifact, which often appears in QSR 
for the exotic hadrons. We show that imposing correct 
criteria on the Borel window rejects such an artifact. In 
Sec lII CI we examine the width effects on the Borel sta- 
bility plots for the effective mass and residue. The way to 
determine the threshold parameter s t h is also discussed 
in detail. 



A. Basic concepts of QSR 

Following the standard way of the QCD sum rule, we 
start with the time-ordered two-point correlation func- 
tion of a tetraquark interpolating field J for the scalar 



n{q 2 ) =i Jd A x e lq - x (0\T[J(x)J f (0)}\0), 



(2) 



where (0| ■ ■ • |0) denotes a vacuum expectation value. 
(Hereafter we write it as (• ■ • ) for brevity.) The QSR 
is then obtained through the dispersion relation, 



poo 

Ren(g 2 ) = P / ds - 
Jo 7T 



1 ImII(s) 



(3) 



satisfying the the spectral conditions, ImJI(s) > 0. For 
sufficiently large — q 2 , the left hand side of can be 
expressed by the operator product expansion (OPE) with 
Ci including the vacuum condensates: 



nr(? 2 ) = E^ (g 2 ) 4 - J iog(-g 2 

3=0 



8+2j 



(4) 



The OPE starts from (g 2 ) 4 log(— q 2 ) reflecting the large 
number of quark fields in the tetraquark interpolating 
field. This turns out to be the main origin of QSR arti- 
facts as discussed in Sec lIIBl 

The imaginary part in the right hand side of Eq.Q 
is deemed to be the hadronic spectrum. Based on the 
quark hadron duality ansatz, we approximate higher en- 
ergy part of the spectral function than a threshold s t h to 
the spectral function obtained by OPE: 

lmU h {s) = 6{s th - s) lmR h {s) + 6(s - s th ) Imn ope (s).(5) 

Here we introduce the threshold parameter sth where the 
quark hadron duality ansatz begins to work. (We use 
the notation E t h = y/s t h in the following.) Hereafter, we 
write the first term of the right hand side of Eq. (O as 



for the later convenience. In usual QSRs, the low energy 
part Il < is parametrized by a delta function, \ 2 6(s — m 2 ) 
with the pole mass m and the overlap residue A of the 
interpolating field J and the hadronic state. Instead, we 
do not specify the form of n< at this stage, since we also 
consider more general contributions from the resonance 
width and background scattering. 

Substituting the expressions Eqs.(U]) and into 
Eq. ([3]) and using Eq. |T]) , we obtain the Borel transformed 
expression of the sum rule 



/ ds e -/M 2 i Imn <( s ) = J2 

JO ■ 1 



(-y c 8+2j - 
j m (M 2 y-i 



poo POO \ 

/ -/ )ds< 

JO J sth / , =0 



s 4 - j , (7) 



with the Gamma function T(n) = (n — 1)!. Calculating 
the Wilson coefficients in the right side, we can evaluate 
the hadronic parameters in the left side as outputs. Here 
the Borel transformation impoves the OPE convergence 
by factor 1/T(j). At the same time, this transformation 
reduces the contaminations from high energy states due 
to the exponential factor e _s / M . 

Using the equality ([7]), the effective mass and residue 
arc derived as 



m 2 eS (M 2 
\ 2 cS (M 2 )^e 



1\J2) = ml a (M 2 )/M 2 



f° th ds e~ s / M2 s Imn<(s) 
f° th ds e-*/ M2 Imn<(s) ' (8) 

tH ds e~ s ' M ' 2 Imn<(s),(9) 



IL < (s) = e(s th -s)Jmn h (s) 



(6) 



in the similar way to mass evaluation in lattice calcula- 
tions. The reason that we call these values as "effective" 
mass and residue is that m c g and A c ff include not only 
pole mass contribution but also the effects of the width 
and the background. In the usual one peak ansatz, the 
width and background effects are regarded as small, and 
then the physical variables should behave independently 
on the Borel mass and Sth is chosen to satisfy this crite- 
rion. (More detailed discussion for how to fix Sth is given 
in ScclHCn 

Now we turn to the explanation for the situation where 
QSR is workable. In any estimation of the physical quan- 
tities, it is required that Eq.© holds with good accu- 
racy. This is satisfied after the suitable selection of the 
Borel mass which achieves the good OPE convergence 
and reduces the unwanted contaminations from high en- 
ergy states. Such a region for the Borel mass is called 
Borel window. This is important both conceptually and 
practically to derive definite results. A conceptual im- 
portance is already emphasized just above. Practically, 
the Borel window is useful to reject the artifacts. An 
explicit example will be shown in Sec lII Bl 

The Borel window in our analysis is determined as fol- 
lows based on Ref. [l(| : The lower boundary of the win- 
dow is set up so as to make the OPE convergence suf- 
ficient in higher-dimensional operators. The criterion is 
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quantified so that the highest-dimensional terms in the 
truncated OPE are less than 10% of its whole OPE, i.e., 



A(M 2 



dim. n terms 



OPE summed up to dim.n 



<0.1. (10) 



At the same time, the higher boundary of the window is 
fixed by the pole-dominance condition that 

f Sth ds e~ s / M2 ImiT < fs) 

Setting up the Borel window is the most important 
step for the application of the sum rules to discuss the low 
energy side in the spectral integral, especially for exotic 
hadrons, as emphasized in [2p| . However, this important 
step was sometimes neglected in the literature for the ex- 
otics. Therefore their successes seem to be partial. Out- 
side of the Borel window, we are stuck with the sum rules 
artifacts, i.e., the artificial stability of the physical quan- 
tities against the variation of the Borel mass. Further, 
the evaluation of the physical properties, mass, residue 
and so on, depends on the selection of Sth so strongly that 
QSR loses the predictive power. In the next subsection, 
we will illustrate the origin of the artifacts using some 
examples. 



B. Pseudo peak artifacts 

In the following arguments, we consider one of the sim- 
plest spectral functions having only the dim.O term pro- 
portional to s 4 in the tetraquark operator case. This case 
study provides a good example to examine typical arti- 
facts in the QSR for the exotics. Similarly this is the 
case if we consider only the lower dimension terms pro- 
portional to s n (n > 0), and then we recognize the im- 
portance of the higher dimension terms dim. 8, 10, 12,... 
beyond the polynomials in OPE. 

Equation (J7|| with only dim.O term becomes 

[ Sth ds e~ s / M2 -Indies) = n'ds e-^'Cos 4 , (12) 

JO T JO 

where the value of the coefficient Co is irrelevant in the 
following discussion. We expect that the dim.O term 
poorly contains the low energy correlation, because the 
polinomial s 4 drastically decreases in small s region, 
while 8 increases rapidly with larger s. Then if we use 
such a simply increasing function as the spectral func- 
tion for the sum rules, we encounter unphysical Borel 
stability for the mass and residue, since the function cut 
above Eth behaves like a peak function as seen in the 
upper left panel in Figfl] We call this pseudo peak art- 
fact which often appears in the QSR without including 
enough higher dimension terms beyond the simple poly- 
naomial s™ (to > 0). 

Let us see how pseudo peak artifact affects the QCD 
sum rules. As seen in the lower panels of Figfl] both 



effective mass and residue plots show the stability in the 
larger M 2 region, as a result that the spectral function 
just below Eth behaves like a peak without large sup- 
pression from the factor e~ s l M . It is important point 
that they exhibit the strong sensitivity to the threshold 
parameter. This is a consequence that the E t h value di- 
rectly determines the position of the pseudo peak. For 
the effective mass illustrated in lower left panel of Figfl] 
the change of E t h by 200 MeV leads the mass change of 
typically 150 MeV. The strong change is unphysical be- 
cause E t h, in principle, does not have any direct relations 
with the position of the resonance peak. The threshold 
dependence is more clearly seen in the effective residue 
plot. Since the spectrum increases like s 4 , taking larger 
E t h leads the drastic increase of the pseudo peak strength, 
as illustrated in lower right panel of Figfl] 

This artifact can be easily found out by the exami- 
nation of the rate of the pole dominance even when we 
treat more complicated spectral functions than that we 
use now. Shown in the upper right panel of Figfl] is the 
pole dominance ratio B(M 2 ; Sth) defined in Eq. tfTTj) as 
the function of M 2 , which is typical quantity to measure 
the strength of the low energy correlation compared to 
the unwanted high energy correlation. The exponential 
factor in Eq. (fT2")l enables to extract the low energy part of 
the spectral function. The ratio B(M 2 ; sth) determines 
the upper bound of the applicable range of the sum rule 
in terms of the Borel mass (Borel window). As seen in 
two lower panels of Figfl] the artificial stability is seen 
above the upper bound of the Borel window, which is in- 
dicated by arrow in the figure. In this way, the condition 
of the sufficient pole dominance rejects the results with 
artificial stabilities. This is one of the practical usages of 
the Borel window. 

Throughout this section, we have seen the origin of the 
QSR artifact outside of the Borel window. The lessons we 
can learn from these examples are the terms with higher 
dimension like dim. 8 have crucial importance to include 
low energy correlation and to be free from artifacts. In 
connection with the criterion on the Borel window, the 
inclusion of higher dimension terms improves the pole 
dominance and considerably extend the upper bound of 
the Borel window. In addition to the ensurance of the 
OPE convergence, this is another reason that we include 
the higher dimension terms up to dim. 12. 



C. Possible width effects on the plot of the 
physical quantites and examination for the threshold 
fixing criterion 

Most of the light scalar nonets are considered to have 
large widths typically 100 ~ 400 MeV though they have 
not been well-determined yet. It is important to check 
how their widths affect the effective mass and residue 
plots. In this section, we employ as a test function a 
simplest spectral function including the width, i.e., Brcit- 
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FIG. 1: (Color online) The examples of the pseudo peak artifacts. The upper left panel is the plot of the spectral function 
with the dim.O term of the OPE up to the thresholds, E t h — 0.8, 1.0 and 1.2 GeV. The rapid enhancement of the polynomial 
s 4 mimics the peak just below s t h = E 2 h . The upper right panel shows the rate of the pole dominance, B(M 2 ; sth)- Without 
higher dimension terms, the rate is very small because of the lack of the low energy correlation. The lower panels show the 
mass in the left panel and residue in the right panel, respectively. The arrows indicate the M 2 values where the pole dominance 
ratio begins to fall below 50%, which is the upper bound of the Borel window. In the larger M 2 , both quantities show the 
moderate stability because of the pseudo peak behavior just below s t h- 
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FIG. 2: (Color online) The relations among the Breit- Wigner form of the spectral function, the effective mass and the residue 
plots. The left panel shows the Breit- Wigner spectral function with the fixed mass 600 MeV and the values of the width are 
taken as 0, 50, 100, 200 and 400 MeV. The middle and right figures show the effective mass and residue calculated by Eqs.(j8]) 
and (J9J with the Breit- Wigner spectral function, respectively. Because of the low energy tails, the effective mass and residue 
plots shift downward in the case of the small M 2 . 



Wigner function (BWF), 



const, x 



r 2 /4 



i^fs-m) 2 +r 2 /4' 



(13) 



where m is the pole mass and V is the width. The const, 
will be taken to normalize the integral of this function to 



1. We show in the first panel of Fig[2]thc BWF with the 
pole mass fixed to 600 MeV and the width taken as 0, 
50, 100, 200, and 400 MeV. 

Let us see the behavior of the physical quantities as a 
consequence after the substitution BWF into Imll < (s) of 
the Eqs.flS]) and ©■ The middle and right panels in Fig [3] 
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show the effective mass and residue plots, respectively. 
In the case of the zero width approximation, both mass 
and residue plots show the complete independence on the 
variation of M 2 as expected. On the other hand, in the 
nonzero width cases, BWF has the tail both in the lower 
and higher energy around the pole mass, and then the 
effective mass shifts downward (uppcrward) in the lower 
(higher) M 2 . 

Fig [2] tells us that, if we consider the effect of the width 
in the effective mass plot, the best Borel stability does 
not necessarily give the best fit of the spectral function. 
This violation of the Borel stability due to the width 
affects how to select the E t h and the final results. In the 
zero width or very narrow width case, usually we can fix 
E th to give the best Borel stability because the effective 
mass plots of the zero width case do not depend on the 
Borel mass M 2 . But in the case with a wider width, how 
should we fix E t h? 

If we assume the shape of the spectrum for the res- 
onances, the effective mass plots obtained by QSR can 
be compared with those of the known spectrum, and we 
can choose E t h to achieve the better fit. This can be 
seen in the following example. Let us consider the p me- 



son which has the mass m r 



770 McV and the width 



r p = 150 MeV, and assume that the spectral function is 
given by the Breit-Wigner form. We show in FigJ3] the 
effective mass plots for the BWF with m p = 770 MeV 
and r p = 0, 150 MeV, and also the effective mass plot 
obtained by the QSR with the OPE given in the litera- 
ture HI] with several thresholds E th = 1.1, 1.2, 1.3 GcV. 
In the zero width approximation, E t h = 1.1 GeV is pref- 
ered within QSR analyses because the Borel stability is 
achieved in the effective mass plot best. In this case, the 
mass is read as 0.7 GeV from the Borel stability. Once 
the effects of the width is considered in the effective mass 
plot, we lose the reason to select E t h which gives the best 
stability. But, comparing the effective mass plot of the 
BWF with the 150 McV width, we find that the case of 
E t h = 1-2 GcV, in which the stability is slightly worse 
than the E t h = 1.1 GeV case, is better to reproduce the 
effective mass plot for BWF with the 770 MeV mass and 
the 150 MeV width. In this way, we can estimate the 
threshold E t h and the resonance mass m. Nevertheless, 
for the scalar meson case, we do not know neither the 
shape of the spectral function nor the width of the res- 
onance, a priori. Thus, we cannot definitely obtain the 
exact values of the threshold and the resonance mass. In 
this analysis, the resonance masses are estimated by as- 
suming the zero width and the 400 MeV width with the 
BWF. Then we adopt the range between two masses as 
our result of QSR. As we will see below, the difference of 
these two mass is about 20%. 

As seen in the above argument, for resonances with a 
wider width, it is more difficult to extract the resonance 
mass than the zero width state. The results unavoid- 
ably have some ambiguities from the assumption of the 
shape of the spectral function and the width. Never- 
theless, we can still estimate physical quantities because 
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FIG. 3: (Color online) The effective mass for the Breit-Wigner 
function with m p =770 MeV, T p =0, 150 MeV, and for the p- 
meson SR results with Eth=l.l, 1.2, 1.3 GeV. 



of the following two points: (i) As illustrated in Figs[2] 
and [21 although the width effects violate the stability of 
physical quantities, the stability is still moderate than 
naively expected even if we take relatively large widths, 
400 MeV. Thus, even if the width exists, the usual deter- 
mination of Eth may work as a first approximation, (ii) 
As an experience from QSR for mesons and baryons with 
small widths, QSR employing E t h taken around the po- 
sition of the second resonances does satisfy the least sen- 
sitibity criterion. Since it is known experimentally that 
the continuous cross section begins from around second 
resonance, to fix E t h around the second resonance is rea- 
sonable and seems to be consistent with the philosophy 
of the quark-hadron duality. Therefore, wc can expect 
the reasonable range of E t h using the experimental facts. 



III. SU(3) SINGLET AND OCTET STATES IN 
THE CHIRAL SU(3) LIMIT 

The a meson may be realized as admixture of the fla- 
vor singlet and octet states, and the /o(980) meson will 
be its nonet partner of the flavor mixing in the flavor 
SU(3) symmetry breaking. It is important for the stud- 
ies of the a and /o(980) meson to know whether both 
resonances should be reproduced in the same footing and 
what makes the differences. For this purpose, we inves- 
tigate quark flavor contents of the correlation functions 
for the flavor singlet and octet states in the SU(3) chi- 
ral limit, in which they decouple each other and we can 
investigate these states separately. In the SU(3) chiral 
limit, we will find that the splitting of the singlet and 
octet states stems only from the annihilation diagrams 
in which some quark lines are not connected between the 
space-time points x and in the intermediate states of the 
correlation function. There are two types of the annihi- 
lation diagrams in the tctraquark case, the color singlet 
two quark state, (qq)i, and the color octet quark with 
gluon, (qq)gG. Wc will emphasize that the contributions 
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from the qq annihilation diagrams are more important in 
the tetraquark correlator than in two quark meson cor- 
relator. 

A. Quark contents of the singlet and octet states 

To study the quark contents of the light scalar nonets, 
it is convenient to introduce the diquark basis, 

U = (ds), D = (su), S = {ud), (14) 

which are anti-symmetrized in the color and flavor spaces. 
We do not specify the Lorentz structure in this section, 
since the flavor contents are the main issue here. The 
details of the interpolating fields will be given in Sec lIV Al 
In the tetraquark case, the quark content of the singlet 
state is described in the diquark basis as 

S= — (UU + DD + SS), (15) 
v3 

and for the octet states, for example, the quark contents 
are given by 

Ox = UD 

e> 2 = -j=(UU + DD- 2SS), (16) 
V6 

in a similar way to the quark contents of the usual qq 
meson cases. 

The isodoublct k and isovector a belong to purely the 
octet because of the nonzero isospin, while the isoscalar 
a and fo can be composed of the mixture of the singlet 
and octet states in the real world where the flavor SU(3) 
symmetry is broken by the quark masses. If the ideal 
mixing is realized, the a and fo can be written as 

o ~ (ud)(ud) = SS=^S~ y|o a , (17) 

fo ~ -^[( us )(us) + (ds)(ds)] 

= -^j=[DD + UU] = \Hs+ sj\02. (18) 

From now on, we consider the chiral limit taking the 
zero quark masses. This limit is convenient to study the 
difference between the singlet and octet states. The study 
of the sigma meson in this work is valid even in the phys- 
ical world with the finite strange quark mass, since, with 
the SU(3) breaking due to the strange quark mass, the 
ideal mixing is expected to be realized and the quark 
contents of the sigma meson is given by Eq. (|17p . In this 
case, the sigma meson consists only of the up and down 
quarks. For the fo and the octet mesons, it is important 
to include the finite strange quark mass. 

We will find that important contributions for the differ- 
ence between the singlet and octet states are the mixing 
of four-quark (4q) states with two-quark (2q) states, two- 
quark and gluon mixed (2qG) states or glueball states. 



B. The annihilation diagrams in 2q-2q and 4q-4q 
correlators 

As seen in Eqs. (fl"7j) and (fl8|) . the difference of the sin- 
glet and octet correlation functions is important to make 
the splitting of the a and fo mesons. Here we would 
like to focus on the annihilation diagrams, since they are 
only the source of the splitting. We will also see that the 
annihilation diagrams play key roles especially for the 
scalar meson spectra. The remarkable point is that the 
effects of the annihilation diagrams are more important 
in the 4-quark correlators than the 2-quark mesonic cor- 
relators. This can be understood from the viewpoint of 
the momentum conservation. Here the annihilation dia- 
gram means the diagrams which have some quark lines 
disconnected between x and 0. 

First of all, we start with discussion about the annihila- 
tion diagrams in the two-quark (2q-2q) correlators. Some 
typical OPE diagrams are shown in FigdJ The dashed 
circles represent a pair of quarks to form the condensate. 
The perturbative diagrams shown in Figf5]a) and b) con- 
tribute relatively high energy side of the correlation func- 
tion. Diagram a) includes 1-loop without a s corrections. 
Diagram b) is one of the annihilation diagrams, since the 
quark lines are not connected in the intermediate state. 
This diagram has three loops with a 2 correction and con- 
sequently it is strongly suppressed compared to diagram 
a) . In diagram b) the two-gluon propagation is necessary 
because of the traceless property of the color SU(3) ma- 
trix. The suppression of this type of the diagrams can be 
also explained from large N c counting [23[. The higher 
dimension terms shown in Figf?] c) and d) contribute to 
the low energy side of the correlation funciton. Diagram 

c) is not the annihilation diagram since the quark lines 
are connected through the quark condensates. Diagram 

d) is very similar with diagram c), but the pairs of the 
quarks to form the condensate are different. Diagram d) 
is the annihilation diagram. However, this contribution 
vanishes at this 0(a s ) order because of the color traceless 
properties. All of these facts suggest the strong suppres- 
sion relative to the non-annihilation diagrams from low 
to high energy in the case of the 2q-2q correlators. 

Next we turn to the discussion of the annihilation dia- 
grams in the four-quark (4q-4q) correlators. We show in 
Figj5]some examples of the diagrams in the 4q-4q correla- 
tors. Diagrams a) and b) are perturbative diagrams. Di- 
agram b) is the annihilation diagram but it is suppressed 
by a 2 in the same way as the 2q-2q correlator case. Di- 
agram c) is the leading annihilation diagram with (qq)i 
configuration. In contrast to the 2q-2q correlator cases, 
the qq annihilation diagrams can appear without gluon 
propagations, since the other two quarks can carry the 
momentum. Consequently no a s suppression comes in. 
Diagram d) is also the annihilation diagram in dim. 8 with 
the configuration (qq)sG. This is not again suppressed 
by the a s correction. Such annihilation diagrams can 
widely appear in the higher dimension terms. This means 
that the annihilation processes are much more important 




a) 





c) 



FIG. 4: (Color online) Some diagrams in the 2-q correlators. 
The flavor is distinguished by the line of color, a) The leading 
perturbative contribution in dim.O term, b) The perturbative 
qq annihilation diagrams in dim.O term with suppression fac- 
tors, 1/N C and a 2 8 . c) The leading power correction in dim. 6 
including (qq)sG configurations without annihilations, d) The 
0(a s ) power correction in dim. 6 to the annihilation diagrams 
is vanished because of traceless properties of Gell-Mann ma- 
trix. 





a) 



b) 





c) 



FIG. 5: (Color online) Some diagrams in the 4-q correlators. 
The flavor is distinguished by the line of color, a) The leading 
perturbative contribution in dim.O term, b) The leading qq 
annihilation diagram in dim.O term, c) The leading power cor- 
rection in dim. 6 to the annihilation diagrams including (qq)i 
configuration without a s suppression factor, d) The anni- 
hilation diagrams in dim. 8 which includes the configuration 
(qq) 8 G. 



in lower energies where we would like to investigate the 
properties of the scalar mesons. 

Here we would like to stress that the annihilation pro- 
cesses do not directly imply 2q propagation emerged from 
tctraquark configuration. Indeed, as shown in Fig[51 the 
annihilation diagrams can represent either 2q s-channel 
propagation or 4q propagation with exchanging the res- 
onance in the t-channel. Especially, contributions from 
the latter processes can be interpreted as diquark-diquark 
correlation in our tctraquark fields as product of diquark 
fields. We will see later that the annihilation diagrams 



annihilation 





s-channel 



t-channel 



FIG. 6: (Color online) The two ways of interpretation of the 
annihilation diagrams. By deforming the quark line, we can 
interpret the annihilation diagrams as either 2q like s-channel 
propagation or 4q like t-channel propagation with exchange 
of mesonic resonances. 



provide imporant contributions in our Borel analyses. 



C. The splitting between singlet and octet states 

Now we discuss how the annihilation diagrams studied 
in the previous subsection contribute to the flavor singlet 
and octet correlation functions. The flavor contents of 
the singlet and octet correlators are given respectively 
by 
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(T[UU(x)] [UU(0) + DD{0) + 55(0)] ) + (perm.), (19) 



(T[UD(x)] [UD(0)]), 



(20) 



where for the octet we use Eq. ([TB| . These correlators 
consist of three types of correlators in terms of the flavor 
content: 



[UU(x)] [UU(0)] 
[UU(x)] [DD(0)] 
[UD(xj\ [UD(0)] 



= [dsds(x)] [dsds(0)] , (21) 
= [dsds(x)] [susu(0)] , (22) 
= [dssu(x)] [dssu(0)] . (23) 



In the first expression, all the diquarks have the same 
flavor, the second expression is the correlation between 
the different flavor tetraquarks consisting a pair of the 
same flavor diquarks, and the third expression has two 
same tetraquarks which are made up of a pair of different 
diquarks. For the flavor singlet, the first two contribute, 
while for the octet only the last type of the correlators 
takes place. 

Now we count possible configurations for the annihi- 
lation diagrams in each correlator shown above. Since 
neither QCD interaction nor vacuum condensates can 
change quark flavors, we can connect only the pair of 
the quark and antiquarks which have the same flavor. 
Thus, the first correlator (f2"Tj) has both the four-quark 
connected diagrams like diagram a) in Fig [5] and the qq 
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annihilation diagrams like diagram c) and d) in Figj5l 

dsds — > dsds (direct), (24) 
dsds [(dd) 1: (dd) 8 G, G 2 } -► dsds, (25) 
dsds -> [(ss)i, (ss) 8 G, G 2 } -> dsds, (26) 

where we explicitly write down the possible color singlet 
configurations for the annihilation processes in the in- 
termediate states. Here we have also mentioned the qqqq 
annihilation process as G 2 in the intermediate state. This 
process gives 4q-glueball correlation. In leading a s anal- 
ysis up to dim. 12, however, such qqqq — > G 2 diagrams 
do not appear. In the similar way, we find that the sec- 
ond type of the correlator (f2"2")) has only the annihilation 
diagrams: 

dsds — ► [(ss)i, (ss)sG, G 2 } — ► susu. (27) 

The last correlator, which contributes to the octet, has 
both the four quark connected diagrams and the qq an- 
nihilation diagrams: 

dssu — > dssu (direct), (28) 
dssu — > [(du)i, (dM)gG] — > dssu, (29) 

As we have mentioned, the first and second types of 
the correlators contribute to the singlet state, while the 
octet correlator has only the third one. From these obser- 
vations, now we can see the singlet state has much more 
annihilation diagrams, or 4q-2q correlation effects than 
the octet states. In fact, the explicit calculation given in 
the AppendixfAl shows the number of the qq annihilation 
diagrams is larger in singlet than in octet by factor 4. 



IV. THE BOREL ANALYSES 

In this section, we perform the Borel analyses for the 
a meson and also the scalar states with the flavor sin- 
glet and octet in the flavor SU(3) limit. We first discuss 
the interpolating fields used in these analyses in subsec- 
tion |IV3 There we introduce linear combinations of the 
interpolating fields constructed by the scalar and pseu- 
doscalar diquarks with a mixing angle 9. We also discuss 
the criterion to determine the interpolating fields based 
on the argument for the Borel window given in subsection 
III Al In subsection IIVB1 we show the results for the fla- 
vor singlet and octet states. Subsection IIV CI is devoted 
to discussion for the result of the a meson. There we will 
find that the correlation function in the Borel analysis 
shows large contributions around E = 600 ~ 800 MeV. 

Throughout our analyses, we use the standard val- 
ues of the condensates in OPE: (qq) = -(0.230 GeV) 3 , 
(qg s aGq)/(qq) = 0.8 GeV 2 , (^G 2 ) = 0.012 GeV 4 and 
a s (lGcV) = 0.4. The OPE is evaluated within the fac- 
torization hypothesis without quantum loop corrections. 
The detail of the OPE calculation is given in AppendixfAl 



A. Interpolating fields 

In QSR, it is technically important to use such good 
interpolating fields as to pick up largely the resonance 
state and have smaller correlation with the higher energy 
states. With the number of quarks increasing, there are 
several choices of the local interpolating fields with the 
quantum number of our interest in contrast to the qq me- 
son cases. The authors of Ref.[22| comprehensively stud- 
ied the possible local interpolating fields of the tetraquark 
and the properties of the scalar mesons by constructing 
QCD sum rules with OPE up to dim. 8. In this work, 
we calculate the OPE of the correlation function up to 
dim. 12, considering linear combinations of the following 
two operators for the a meson: 

Jp = e abc e dec [ulCd b }[u d Cd^}, (30) 
Js = e abc e dec [u^Cj 5 d b ][u dl5 Cd^}, (31) 

where a,b,c... represent the color indices, i\) T means the 
transposed spinor of ip an d C is the charge conjugation 
matrix. The above interpolating fields, Jp and Js, are 
constructed by two pscudoscalar diquarks and two scalar 
diquarks, respectively. The linear combination is given 
with a mixing angle 9 (0 < 9 < tt) by 

J (9) = cos 9 J P + sin 9 J s . (32) 

The mixing angle 9 in the interpolating field (f3"2"| is 
determined so that the interpolating field couples more 
strongly to the resonance state and have less contribu- 
tions from the higher energy states and scattering back- 
ground, which is achieved by imposing the following cri- 
teria: 

a) . Sufficiently wide Borel windows are established in 

the sum rule. 

b) . The effective masses and residua are weakly depen- 

dent on the Borel mass M 2 as much as possible. 

c) . The results have also adequately weak dependence 

on the threshold E t h- 

d) . The effective residua are satisfactorily large. 

These criteria are not independent each other and the 
mixing angle can be insensitive to some of them. Hence, 
the mixing angle is not uniquely fixed in quantitative 
manner. Here we discuss typical mixing angles which 
simultaneously satisfy the above criteria with better de- 
gree. In the followings, we explain the meanings of the 
criteria separately and understand the priorities among 
these criteria. 

The criterion a) is most important, being essential to 
avoid the pseudo peak artifacts and reduce truncation 
errors of the OPE, as emphasized in Sec lII Al Thus, this 
criterion should be always satisfied independently of the 
other criteria to guarantee the sum rules contain largely 
contributions from low energy physics. 
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After we find the Borcl window in the sum rules, the 
next task is to establish better pole isolation from the 
background contamination, which can be done by satis- 
fying the criteria b) and c). The background contamina- 
tion comes mainly from scattering states of two mesons 
and other resonance states. In the low energy side, scat- 
tering states of two Nambu-Goldstone bosons, i.e., irir, 
KK, rjn, scattering states can appear from E = in 
the present analyses because all the Nambu-Goldstone 
bosons are massless in the chiral limit. If these states 
largely contribute to the sum rules, the effective residue 
and mass obtained in the sum rules have larger Af 2 de- 
pendence in the M 2 region below the resonance pole than 
we expect from the resonance width, which has been dis- 
cussed in Sec lII CI This contamination below the res- 
onance pole can be reduced by the selection of 9 with 
satisfying the criterion b) well. The background contam- 
ination above the pole energies comes from higher excited 
resonance states and meson scattering states other than 
the Nambu-Goldstone bosons. In the case that such con- 
tamination is large, the physical quantities have large 
dependence against the variation of Eth- This can be 
reduced by imposing the criterion c) . 

The criterion d) is imposed to have the large overlap 
of the interpolating field with the resonance pole and to 
reduce the truncation errors of OPE. In our analyses, 
we neglect higher-dimensional terms of OPE than the 
dim. 12 terms. If the residua obtained in the sum rule 
arc sufficiently large, the interpolating fields used in the 
analyses have enough overlaps to the resonance states 
and the sum rules constructed with the approximated 
OPE may contain so much information for the resonance 
to be extracted. But, if the obtained residua are small, 
the neglected higher-dimensional terms might have im- 
portant information on the resonance states, although 
the OPE convergence is well provided. Thus, we con- 
sider the criterion d). Nevertheless, the criterion d) has 
less priority than the criteria b) and c) . This is because 
the best choice of 9 is not always the one which have the 
largest overlap with the region below E t h- Even if the 
low energy correlation is large, it can be just a signal of 
strong contamination from the background. Therefore, 
the most important thing is that we construct the sum 
rules within the sufficiently wide Borel window satisfying 
the criterion a) with the best choice of the interpolating 
field which have the least overlap with the background. 

B. Borel analyses of flavor singlet and octet states 

Let us consider the flavor singlet and octet states pro- 
duced by the interpolating fields given in Eq. (|32[) with 
the mixing angle 9 and having the quark contents (fT5|) 
and (|16|) . The mixing angle is determined by the criteria 
discussed in Sec lIV Al The detailed discussions and the 
technical issues are given in Appendix|5J Our findings 
for the mixing angle are as follows: 

• the Borel windows can be established in all the mix- 



ing angles except 9 ~ for the flavor octet case. 

• the results for the physical quantities are not very 
sensitive to the choice of the mixing angle 9 except 
some region of 9 both in the flavor singlet and octet 
cases. 

• one of the best value of the mixing angle both for 
the singlet and octet is 9 ~ 77r/8. Thus in this 
section we will show the results for 9 = 7ir/8. 

• the residua of the singlet states have larger values 
than those of the octet states in almost all region 
of 9. 

• the effective masses of the octet case are typically 
smaller than those of the singlet, and the Borel sta- 
bility is worse in the octet case. 

The last two findings may imply, in comparison to the 
singlet case, that the octet interpolating field used here 
does not have sufficiently large overlap with the reso- 
nance state or that the width of the resonance in the 
octet channel is considerably large. For further investi- 
gation of the octet states, it might be necessary to use 
interpolating fields with different Lorenz structure. 

Hereafter we proceed to discussions of the Borcl anal- 
yses of the scalar mesons with fixing the mixing angle as 
9 = 7it/8. First of all, we examine the OPE convergence 
of the pole dominance, which are important to establish 
the appropriate Borel window, and we emphasize impor- 
tance of the higher dimension terms. Here we fix the 
threshold to a typical value E t t = 1.0 GeV. We show in 
FigEJthc OPE dependence of the ratio of the pole contri- 
bution to the whole spectral function B(M 2 ; sth) defined 
in Eq. pTjl for singlet and octet cases. We can see that 
the OPE corrections to the perturbative term (dim.O) is 
extremely important to achieve the pole dominance. It 
can be also seen in Fig[7]that the higher dimension terms 
than dim. 8 is not so important for the pole dominance. 
Nevertheless, as seen in the plots of the effective masses 
(Fig|5]), the dim. 10 and 12 terms are important to ob- 
tain the Borcl stability in the small M 2 region (typically 
A/ 2 < 1.0 GeV 2 ~ A^m ax ) where we will evaluate the 
physical quantities later. This means that these higher 
dimension terms have essential to reproduce the low en- 
ergy resonances. For the singlet case, the dim. 12 term 
improves the Borel stability well, while in the effective 
mass for the octet state the dim. 10 and 12 terms do 
not help as well as the singlet case. The Borcl stabil- 
ity obtained with including only the dim.O term can be 
interpreted as the artificial stability discussed in Scc lIIBl 
since the pole dominance is not achieved in these cases. 
As emphasized in Sec lIII Cl the difference of the effective 
masses of the singlet and octet comes from the annihila- 
tion diagrams, which is increasingly important in higher 
dimension terms or low energy region. 

Next we evaluate the resonance masses and residua 
for the flavor singlet and octet states in the chiral limit. 
We calculate the masses and residua widely with several 
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FIG. 7: (Color online) The OPE truncation dependence of 
the pole dominance, B(M 2 ; sth), for singlet and octet case 
with the threshold E th = 1.0 GeV. 
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FIG. 9: (Color online) The effective masses for the singlet and 
octet states. The values of the threshold Eth are 0.8, 1.0, 1.2 
and 1.4 GeV. The downward and upward arrows indicate the 
lower and upper bounds of the Borel window, respectively. 
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FIG. 8: (Color online) The OPE truncation dependence of the 
effective masses for singlet and octet case with the threshold 
Eth = 1.0 GeV. 



values of Eth = 0.8, 1.0, 1.2, and 1.4 GeV, in order to 
see the threshold dependence discussed in Sec II V Al and 
the effect of the resonance width studied in Sec lII CI We 
show in Fig|5] the plots of the effective mass for singlet 
and octet states with E th = 0.8, 1.0, 1.2, and 1.4 GeV. 
The downward and upward arrows indicate the lower and 
upper bounds of the Borel window, respectively. The 
way to fix the Borel window has been given in Eqs. (|10p 
and (fTT|) . The lower Borel mass is not dependent on 
the threshold, since it is determined only by the OPE 
convergence. 

We extract the masses of the scalar mesons from the 
effective mass plot in FigOwith consideration of its possi- 
ble modifications from the resonance width, as discussed 
in Sec lII CI For the singlet state, we find the moder- 
ate Borel stability for E th = 1.0 ~ 1.4 GeV in FigH 
This would imply that some resonance state saturate the 
spectral function around E ~ 700 McV. We estimate 
the mass for the singlet state at 700 MeV for a small 
width state, finding almost perfect Borel stability with 
Eth = 1-0 GeV. If the state has some larger width, we 
evaluate the mass by 850 MeV with a possible width 400 
MeV according to the discussion in Sec lII CI in which we 
found that the Borel stability is not perfectly achieved in 
the case of the resonance with width. Finally we conclude 
that the present QCD sum rule estimate the flavor sin- 
glet scalar meson by 700 ~ 850 McV with consideration 



of the possible width up to 400 MeV. 

For octet state, we estimate the resonance mass at 
600 ~ 750 MeV in the similar way to the singlet case. 
Nevertheless, the effective mass in the octet channel de- 
pends on the Borel mass fairly. One of the possible ex- 
planations is that the signal of the octet resonance is 
weak and considerably affected by low energy scattering 
states. Another possibility is that the observed state has 
the large width, let us say, 100 ~ 400 MeV, which can be 
estimated from the comparison to the effective mass of 
the Breit-Wigner form in FigO Probablly both effects 
may play important roles in the octet sum rules. 

We also plot the effective residua for the singlet and 
octet cases in FigJlOl The Borel stability is achieved 
as similarly as that of the effective mass. The residua 
for single and octet states are evaluated by (20 ~ 35) x 
10" 7 GeV 8 and (7.5 - 13) x 10" 7 GeV 8 , respectively. 
The overlap strength of the singlet case is factor 2^3 
larger than the octet state. 

Both in the singlet and octet states, the threshold de- 
pendence is not very small. This would suggest the sep- 
aration between resonance pole and threshold would not 
be completely done. For the further investigation, we 
would need more elaborated technique to isolate the reso- 
nance states, for instance, large N c argument as discussed 
in Ref.[I3. 



C. Borel analysis for sigma meson 

We discuss the Borel analysis for the a meson in the 
chiral limit. With the SU(3) breaking, the ideal mixing of 
the flavor singlet and octet components may be realized 
and the lighter scalar meson with 1 = and S = 0, 
namely the a meson, may be described by only the up and 
down quarks. (The ideal mixing is the assumption in this 
work. To check this assumption, we need to investigate 
the flavor mixing angle in QCD sum rules with a finite 
strange quark mass.) For the ideal mixing, the spectral 
function for the a meson is given by a linear combination 
of the singlet and octet spectral functions in the chiral 
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FIG. 10: (Color online) The effective residua for the sin- 
glet and octet states. The values of the threshold E t h are 
0.8, 1.0, 1.2 and 1.4 GeV. The downward and upward arrows 
indicate the lower and upper bounds of the Borel window, 
respectively. 
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FIG. 11: (Color online) The effective mass for the a meson. 
The values of the threshold E th are 0.8, 1.0, 1.2 and 1.4 
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and upper bounds of the Borel window, respectively. 



limit: 

{TJj a ) = ^{TJ s J s ) + ^{TJ Jo). (33) 

We work again with 9 = 7tt/8 for the mixing angle of the 
interpolating field. 

Evaluating the effective mass and residue of the sigma 
meson, we obtain the plots with E^ = 0.8, 1.0, 1.2 and 
1.4 GeV in Figs fTTI and [T2l In these plots, we find again 
fairly good Borel stability both for the effective mass and 
residue. In the same way as the flavor singlet case, we 
estimate the resonance pole mass at 600 ~ 800 MeV with 
paying attention to the possible width ~ 400 MeV. The 
pole residue is also evaluated by (10 ~ 20) x 10~ 7 GeV 8 . 
The value of the residue is consistent with the expec- 
tation that it is a superposition of the singlet and octet 
states with group theoretical weights for the ideal mixing, 
which are 1/3 and 2/3 for the singlet and octet contribu- 
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tions, respectively. 

This a mass obtained in the present sum rule analy- 
sis with the tetraquark operator is considerably smaller 
than that extracted from correlator analyses with the 
two-quark operators, which is ~ 1.2 GeV. As discussed 
in Sec lIV Bl the 2q operator mainly couples to the 2q 
components, while the tetraquark operators induce both 
the 4q components and the considerable 2q and 2qG com- 
ponents. The present tetraquark investigation shows that 
that inclusion of 4q and 2qG components is relevant to 
explain the smaller a mass 24 1 . 



As a whole, finding the wider Borel window in the 
sigma meson sum rule, we conclude that the results for 
the a meson is more reliable than the singlet and octet 
cases. Especially it is notable that the lower bound of 
the Borel mass is sufficiently small. This helps us to 
investigate low energy contributions for the a channel. 
Due to the small lower bound of the Borel mass and the 
Borel stability obtained even in M 2 = 0.3 MeV 2 , we con- 
clude that the pole dominance of the a resonance is fairly 
established with small contamination of low energy scat- 
tering. We also see that the large strength comes mainly 
from the flavor singlet component, which shows the large 
enhancement around E = 700 ~ 850 MeV region. 

In closing this section, we would like to emphasize that, 
although having some ambiguities coming from the influ- 
ence of the scalar meson width on the Borel stability and 
the choice of the mixing angle 9, the present QSR analy- 
ses with the tetraquark operators provide the scalar me- 
son masses less than 1 GeV, which cannot be produced 
in QSR with 2q operators. 
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V. DISCUSSION 

The present QSRs are constructed in the SU(3) chiral 
limit. It is very interesting to extend our QSRs beyond 
the chiral limit with a realistic strange quark mass. Es- 
pecially for realistic calculation of the masses of the octet 
scalar mesons, such as k and ao, inclusion of the strange 
quark mass is absolutely necessary. It is also interesting 
to investigate the mixing angle of the flavor singlet and 
octet interpolating fields with finite strange quark mass. 
In this work, we have assumed the ideal mixing for the 
flavor content of the a meson operator. But it is not 
trivial to realize the ideal mixing in the physical scalar 
meson nonet. Analyses with the strange quark mass for 
the flavor mixing angle is possible within the QSR ap- 
proach. This would clarify the strange contents of the a 
meson. 

Inclusion of the quark masses with the SU(3) symmet- 
ric way is also a good theoretical issue. In the present 
analyses we have found that contributions from the anni- 
hilation diagrams are responsible for the splitting of the 
singlet and octet states. It is good to know the quark 
mass effects on the splitting. Especially, in the chiral 
limit, we have found that the octet tetraquark operator 
weakly couples to the physical state, so that it is inter- 
esting to investigate how the octet state changes with the 
finite quark mass. 

The issue of the flavor mixing angle and splitting of the 
singlet and octet states is closely related to the mass dif- 
ference of the /o and ao- In the present experiments they 
almost degenerate. The ao meson may be classified into 
the flavor octet, while the /o meson may be described by 
a linear combination of the singlet and octet operators 
with the flavor SU(3) breaking. Thus the mass differ- 
ence is given by the mixing angle and contributions of 
the annihilation diagrams in the QSR approach. It is 
interesting to study how the degeneracy of /o and ao is 
realized in QSR. 

In this work, we have used the tetraquark operators 
with expectation to have strong coupling to the four- 
quark states and some contributions from the two-quark 
(2q) and two-quark-one-gluon (2qG) states. It is also in- 
teresting to investigate the scalar meson states using the 
operators having large overlap with two-quark states and 
small coupling with four-quark states and to compare this 
with the present analyses. It has been reported in Rcf. 8] 
that k and ao may have large components of two-quark 
states rather than four-quark states and the coupling of 
the two-quark states with nonresonant scattering states 
makes the k and ao masses lighter in a similar way to the 
Feshbach resonances. In fact, the previous QSR studies 
with two-quark mesonic currents for k(800) and ao(980) 
showed 1.0 ~ 1.2 GeV and this is not too higher than 
the experimental value. It would be natural to develope 
the ideas that the 2q state with some additional compo- 
nents might lead the experimental mass. In future, we 
will report the analysis using such operator with the help 
of the large N c argument [27[ ■ 



In order to investigate further the widths of the scalar 
meson states, it is necessary to investigate three point 
correlation functions for the a — > irir, for instance. To 
obtain the decay width in QSR, we need to calculate both 
the two-point and three-point correlation functions, and 
to perform combined analyses of them in self-consistent 
ways to determine the mass, the residue, the threshold 
and the mixing angle. 

It is also interesting to calculate two photon decays of 
the scalar mesons in QSR in order to compare mesonic 
molecule pictures. It has been reported in Ref.[26[ that 
the KK molecule picture of the /o and ao mesons pro- 
posed by Ref.[25[ provides a factor 3 larger decay rate 
of ao — > 77 than experiments and also that two quark 
picture predicts much larger decay rates of ao — * 77 and 
fo ~ * 77- The discrepancy could be explained by the 
tetraquark picture, since decay of four quarks into two 
photons may be suppressed compared to two quarks to 
two photons by the electromagnetic constant a and the 
quark wavefunction |^(0)| 4 at the origin. 

It is an important question whether correlation func- 
tion analyses with local operators, such as the present 
QSR, can reproduce nature of spatially extended objects 
like mesonic molecular states. The correlation functions 
observe overlaps of the local operators and the physical 
states. In principle, if the overlap is not extremely small, 
the correlation functions have components of the molecu- 
lar states. However, it is not sufficient at all to investigate 
the nature of the physical states only from the overlap of 
the local operators and the physical states. We certainly 
need other information of the wavefunctions, such as the 
decay widths. 



VI. SUMMARY 

In this paper, we discuss the properties of the er(600) 
meson in the QCD sum rules (QSR) for the tetraquark 
operators, emphasizing the importance of setting up the 
Borel window to avoid artificial Borel stability. The fla- 
vor structure of the sigma interpolating field is given by 
the SU(3) singlet and octet representations assuming the 
ideal mixing with the SU(3) flavor breaking. We invesi- 
gate also the flavor SU(3) singlet and octet states for 
the light scalar meson with 1 = and S = in the 
massless limit with special attentions for the roles of the 
constituents of the a meson. Having shown that the split- 
ting between the singlet and octet states stems from the 
annihilation diagrams in which only two quark lines out 
of four carry the hard momenta, we have found that the 
contributions from the annihilation diagrams are respon- 
sible for the mass of the a state around 600 ~ 800 MeV, 
which is much lighter than the mass extracted by QSR 
with the two quark operators. 

Investigation of the scalar meson in the present QSR 
have been carried out using the correlation functions for 
a linear combination of the tetraquark operators made 
up by the scalar and pseudoscalar diquarks. The mixing 
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angle has been determined with QSR so as to achieve 
better resonance isolation from the background contam- 
ination. We have also discussed possible influence of the 
large width of the scalar meson on the Borel stability. 

For the singlet state, having found fairly good Borel 
stability within the practical Borel window, we have ob- 
tained the mass as 700 ~ 850 McV with considerations 
for the ambiguities coming from the possible width of the 
scalar meson and small dependence on the threshold pa- 
rameter. For the octet state, we have found rather poor 
Borel stability and obtained a bit lighter mass around 
600 ~ 750 MeV with smaller coupling strength than 
the singlet case. This would suggest that, in the octet 
channel, the signal of the resonance state is rather small 
compared to background scattering states of two Nambu- 
Goldstone bosons, or that the resonance state has the 
large width. The present analyses have been done in 
the SU(3) chiral limit. It is interesting to investigate the 
octet state beyond the chiral limit with a realistic strange 
quark mass for the the k and clq mesons. 

For the a meson, having found wider Borel window 
and more Borel stability than the flavor singlet and octet 
cases, we have estimated the resonance mass at 600 ~ 
800 McV and the residue at (10 - 20) x 10" 7 GeV 8 , 
which include again the ambiguities of the width and the 
choice of the threshold. Although the present analyses 
have some ambiguities coming from the influence of the 
physical width of the scalar meson on the Borel stability 
and the choices of the threshold and the mixing angle, 
we conclude that the QSR with the tetraquark operators 
provides the resonance masses of 600 ~ 800 MeV in the 
a meson channel in the chiral limit. This suggests that 
the a meson has largely the four-quark components, since 
the QSR with the 2q operators cannot produce masses 
less than 1 GeV in the scalar channel. 
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APPENDIX A: THE OPE RESULTS 

In this appendix, we present the expression of each 
term of the OPE for the correlation function ([2]) with 
the interpolating field (j3"2")l having the mixing angle 9 
(0 < 9 < tt): 

J (9) = cos 9J P + sin 6 J s (Al) 

The interpolating fields with specific mixing angles have 
definite chiral representations: J(3ir/4) is classified into 



the irreducible representation [(3, 3)<g>(3, 3)] of the chiral 
SU(3)l<8)SU(3)h group, while J(n/4) is a combination of 
the [(8, 1) £g> (1, 8)] and (1, 1) representations. (Actually 
J(7r/4) is classified into the nonet [(9, 1) <g> (1, 9)] repre- 
sentation of the U(3)i®U(3)fl; group.) 

The OPE of the correlation function is expressed by 
the the coefficients C,: 

4 oo „ 

nfV) = E (9 2 ) 4 - J lo s(-? 2 ) + E tW- ( A2 ) 

o=o j=i [q ' 

The relevant diagrams to calculate the OPE are shown 
in FigQJ 

We use c = cos#, s = sin6>, and (uu) = {dd) = (ss) = 
(qq) for the simplication of the OPE expressions. We 
also define A as a number of the annihilation diagrams, 
N = 1 (for octet), N = 2 (for a) and N = 4 (for singlet). 



15 



For the terms from dim.O to dim. 4, we have 

1 



C = -(c 
C 4 = -(c 



2 1 s 2 ) 



2 1 s 2 ) 



2 12 3 5tt 6 ' 



2 10 3 TT 4 



x 2, 



(A3) 
(A4) 



where "x..." is the multiplication factor from the permu- 
tation of diagrams. 

From dim. 6, the annihilation diagrams begin to ap- 
pear. The dim. 6 terms are calculated as 



C 6a = (c 2 - s 2 ) 



2 3 3 7T 2 



C, 



6b 



-(c 2 +s 2 - 2sc) 



x2, 
(gg) 2 

2 5 3 TT 2 



x N, 



(A5) 
(A6) 



where the suffix of C stands for the diagram in FigfT31 

The dim. 8 terms provide most important contribu- 
tions. These are the first contributions beyond the poly- 
nomial s n (n > 0) and give the most of the low energy 
correlation. They are calculated as 



C 8a = -(C 2 



(qq)(qg s aGq) 



c 8 



■s~) - x 4, 

3 7T 2 

(qq) (qg s aGq) 



-(c 2 + s - 2sc) 



2 6 3 7T 2 



(A7) 
x 2N. (A8) 



Note that the annihilation diagrams can give substantial 
contributions roughly ~ 40% of the total of dim.8. 

From dim. 10, the OPE begins to converge, but still 
substantial contributions give to the correlation function. 
The coefficients of the dim. 10 terms are obtained as 



C 1Qa = (c 
C\ob = 
C Wc = (c 

cw = 

ClOe = 



2 , (qq) 2 (^G 2 ) 
2 g 2 VM \ ir / x 4 



c 2 -s 



2 ^ 



2 3 3 3 
2 ^ (qg s (TGq) 2 

2 6 3 TT 2 
(qq) 2 (^G 2 



2 4 3 2 TT 2 



x2, 
x 2, 



^ c 2 _ s 2\ n (qgsvGq) 2 



2 10 3 2 TT 2 



x2, 



(A9) 
(A10) 
(All) 
(A12) 



(c 2 + s 2 - 2sc) {qq)2 Jl 2 ° 2) x N, (A13) 

C 10/ =-(c 2 +s 2 -2sc)%^x7v. (A14) 

In our analysis, we calculate the OPE up to dim. 12, be- 
cause further dimension terms have no large enhance- 
ment factor of ~ (47r) 2 coming from cutting quark loops 
[2(il ]. There are many diagrams for the dim. 12 terms be- 
cause we have many ways to attach the gluon line to 
the quark lines. The results after adding the all possible 
diagrams are 



Cl2a 



m 4 



3 2 

x [8(c 2 + s 2 ) - 3N(c 2 + s 2 - 2sc)], 
2 4 7ra s (gV 4 



(A15) 
(A16) 



C 



12b 



3 s 



[2(c 2 -s 2 )-7V(c 2 -sc)].(A17) 



There are no terms having the structure other than the 
above, such as {qq){qg s aGq){ 9 z-G 2 ) . 



APPENDIX B: 9 DEPENDENCE OF THE 
EFFECTIVE RESIDUE AND MASS 

In this Appendix, we discuss the detail of the 9 de- 
pendence. We would like to draw a rough sketch of the 
allowed region for the mixing angle to achieve wider Borcl 
window and better pole isolation following to the discus- 
sion in Sec IIV Al 

We show in FigQj] the effective mass obtained 
by the QSR as functions of the mixing angle 9 
(0, 2tt/8, 5tt/8, 7tt/8) and the Borel mass M 2 . The 
threshold is taken as E th = 0.8, 1.0, 1.2, 1.4 GeV in the 
same way as subsection IIVBI and IIV CI We also show 
the effective residue in Fig[l5]with the same 9 and E t h 
as the effective mass plots. 

The essential points were already mentioned in subsec- 
tion |IVB] Here we give more detailed explanations. 

• the Borel windows can be established in all the mix- 
ing angles except 9 ~ for the flavor octet case. 
For 9 ~ 0, the spectral condition ImlI ope (s) > 
is slightly violated unless we take Eth as large as 
~ 2.0 GeV. Thus the region around 9 = was 
avoided in subsection IIVBI to compare singlet and 
octet cases in the region where sum rules for both 
cases well satisfy the criteria in subsection IIV Al . 

• the results for the physical quantities are not very 
sensitive to the choice of the mixing angle 9 in the 
regions of 9 = ~ 27r/8 and 6n/8 ~ tt both in 
the flavor singlet and octet cases. This is most 
clearly reflected in the effective residue plots in 
FigdSl On the other hand, except these regions, 
the effective residua strongly depend on the mixing 
angle and the threshold. Wc interpret that these 
strong dependences are mainly related to contribu- 
tions from scattering states. This is because the 
residua in these 9 regions show the large E t h de- 
pendence, which indicates that the resonance pole 
is not isolated enough. 

• one of the best value of the mixing angle both for 
the singlet and octet is 9 ~ 77r/8. If we consider 
only the singlet and a meson cases, the best result 
is obtained with 9 = 0. In this case, the threshold 
dependence is extremely small, and we can expect 
that there arc small contaminations in the region 
from resonance pole to E t h, i-c, the resonance is 
isolated. 

• the residua of the singlet states have larger values 
than those of the octet states in almost all region 
of 0. 

• the effective masses of the octet, case arc typically 
smaller than those of the singlet, and the Borcl sta- 
bility is worse in the octet case. 
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FIG. 14: (Color online) The mixing angle (9 — 0, 2n/8, 5n/8, 7ir/8) dependence of the effective mass plots for singlet, octet, 
and a meson cases. For 9 — in the octet case, we show only the lower bound of the Borel window since the spectral condition 
is not satisfied in the region Eth = 0.8 ~ 1.4 GeV. 



In addition to these, we put one remark on the chiral 
representation of the interpolating fields: in the case of 
J(n/4), i.e., the combination of [(8, 1)® (1, 8)] and (1, 1) 
representations, the contributions from annihilation dia- 
grams vanish (see also OPE expression), and the results 
of singlet, octet, and a meson degenerate. 

Taking into account the above all remarks, we evalu- 
ate the mass for singlet, octet, and a meson cases. We 
mainly consider the cases of 6 — and 7tt/8 in FigJI31 



For the singlet case, the typical effective mass ranges are 
from 600 MeV to 850 MeV in wide Borel windows with 
reasonable Borel stability. The E t h dependence is also 
reasonalbly small. For the octet case, we must avoid 
6 ~ 0. The 9 ~ 7tt/8 case was already discussed in 
Sec lIV Bl Finally, for a meson case, the behavior of ef- 
fective mass plot is better than both singlet and octet 
cases. We evaluate the mass as 600 ~ 800 MeV. 
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